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Abstract 

In this article we study the formation of the two modes Bose-Einstein condensate and the cor- 
relation between them. We show that beyond the mean field approximation the dissociation of a 
molecular condensate due to the chemical potential curve crossing leads to the formation of two 
modes condensate. We also show that these two resulting modes are correlated in a two mode 
squeezed state. 

PACS numbers: 03.75.Fi, 03.65.Nk, 34.50.Lf, 34.50.RK 
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INTRODUCTION 



Recent studies on Bose-Einstein condensate have directed the attention of many theorists 
and experimentalists to the multimodes Bose-Einstein condensate. In particular the appear- 
ance of multimodes condensate from a molecular one via molecular condensate dissociation 
process. In |IJ], the authors studied the dissociation of a molecular condensate into a 
one mode atomic condensate. Moreover, Yurovsky et.al showed that the resulting atomic 
condensate is in a squeezed state. They also showed that the number of atoms and 
squeezing parameters in the case of time dependent chemical potentials depend on the rate 
of curve crossing and the coupling between the molecular and the atomic fields. 
The problems of curve crossing can be divided into two different classes. These are the 
r-dependent crossings and the t-dependent ones. In the former the potentials are coordi- 
nates dependent and the crossing takes place in space, while in the latter the potentials 
are time dependent. The standard approach to the crossing problems is the semiclassical 
approach which is known as Landau-Zener (LZ) theory. The LZ theory provide a very 
good description of the crossing problem for systems with potentials that cross at one point 
only. However, if crossing between two potentials occur at two points or more then the 
solution of the problem is possible by using a different approach which is out of the LZ theory. 

In the production of an atomic condensate from a molecuar one the curve crossing problem 
arise from the crossing of the chemical potentials of the condensates. In this article, we 
study the creation of multi modes atomic condensate from a molecular one near the curve 
crossing of the system's chemical potentials. We show that within reasonable approximations 
the problem can be reduced to the case of two atomic modes and one molecular mode 
condensates. We show that exact results for this problem are still possible since the equations 
for the operators evolution decouple so that the two level Landau-Zener problem can be 
exactly solved. Compared to one atomic mode process, the two modes system provides much 
more possibilities for correlated states since the initial conditions depend on two complex 
parameters which describe the initial states of the two modes. The aim of this article is to 
generalize the results of |lj to the case that includes exited modes of the atomic condensate. 
In this article we study the process when the dissociation of atomic collisions in one mode 
leads to the production of atoms in two different modes. 
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THE PROBLEM 



In it was shown that the description of molecular and atomic condensate by two 
modes only is generally not realistic, since the exited atomic modes interact with molecular 
condensate. Consequently some of them become highly populated and must be taken into 
account. Generally speaking the Hamiltonian that includes more than one atomic and 
molecular modes can be written as follows 0: 

H = ^2(fi k ala k + i/fc-ViVk) + 9k,k' ^k+k'^k' + ipk+k'aWw 
k 1 k,k' 

+ -Qint a k+k'-k" a k" a k' a k (1) 

1 k,k',k" 

If we assume that the chemical potentials Hk depend on time then the Landau-Zener model 
for multimode case becomes very complicated to solve even in the case when the molecular 
field is treated as a c-number. However, if the number of atoms in the atomic modes is 
not very large we can disregard the interaction term (g int = 0). Let us consider a system 
that initially has one highly populated mode of the molecular condensate. For example, 
molecular mode with total momentum equals zero. Due to momentum conservation, the 
atomic mode with momentum k would couple via the interaction with the molecular mode 
to the atomic mode of momentum —k only. In such systems the problem is reduced to the 
interaction of two atomic modes with molecular condensate. Other molecular modes can 
also be considered but in this article we will show that the effective coupling constant in 
this case is propotional to the square root of the number of particles in the molecular mode; 
therefore, we can disregard the interactions with the initially empty molecular modes. 

Let us consider a system of two atomic fields A and B and one molecular field AB. 
Suppose a, b and ip are the annihilation operators of A, B and AB field respectively. If we 
disregard the interaction between A and B during curve crossing then the Hamiltonian for 
molecular dissociation in the process AB — > A + B is time dependent and can be written as 
follows: 

H = /^(t)a + a + fx 2 (t)b + b + gip + ab + g*4)a + b + (2) 

where the molecular field energy is set to zero. This Hamiltonian different from the one 
considered in [[J where A = B. Here a and b are anihilation operators of distinct atomic 
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modes so [a + ,5] = 0. Only atomic mode with zero momentum can be described by the 
process A 2 —>■ A + A. In the case of Hamiltonian (|2|) the atomic modes A and B become 
correlated. We assume that there is a macroscopic number of molecules that does not change 
considarably during the process so that we can substitute < ip > instead of ip. Substituting 
g < ip >*= 7 into (0) we get the following atomic Hamiltonian: 

H = Hi{t)a + a + /i 2 (t)b + b + ^ab + 7*a+&+ (3) 
Schrodinger's equation leads to the following operator equations: 

ia = H\{t)a + 7*6 + 
ib+ = -ii 2 (t)b+ - 7a (4) 



and 

ia + = —fii(t)a + — 76 
ib = /2 2 (t)b + Ya + 
If we make the following change of variable: 



(5) 



a _» ae' { S a*i(t)/2-w(*)/2)«K 

I _> le 1 f (Ml(*)/2-M2(*)/2)dt ^ 

f jL(t) = fji 1 (t)/2 + fi 2 (t)/2 



then equations (§),(n) read: 



ia = /i(t)a + 7*6 + 
ib + = —ji(t)b + — 70, 



(7) 



and 



(8) 



ia + = — n{t)a + — 76 
ib = fi(t)b + 7*a + 

In the adiabatic approximation there are no transitions except for the case when adia- 
baticity is violated near curve crossing points fi(t) = 0. According to Landau-Zener theory, 
to estimate transition amplitude we can approximate fi(t) linearly near the crossing points, 
hence the equation that should be considered are: 

ia = Bta + 7*6+ 

' (9) 
ib+ = -ptb + - 7a 
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The corresponding second order differential equation for operator a reads: 

l + ((3H 2 - | 7 | 2 + ^ 2 )a = (10) 

which is exactly the same operator equation that appears in the case of one mode condensate. 
The details of its solution can be found in [l|] . The difference from the one mode process is 
that the operator a in @ is coupled to operator b + rather than a + . The solution in the two 
mode case can be written in the form: 

a(t) = (t) c (t)a(t Q ) + (j) s (t)b + (t ) (11) 

with initial conditions: 

c (to) = l,0.(to) = O 

Mt ) =pt ,4> a (t ) = -17* 

The average numbers of the atoms in the A-mode are: 



< a + (t)a(t) >= n st + n 



sp 



(13) 



where 



n 



sp 



(14) 



corresponds to the spontaneous transitions into the A— atomic vacuum state. This term 
does not appear if the problem is treated in the mean field approximation and is the result 
of quantum effects [[ 
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The quantity n st corresponds to the stimulated transitions in the case when atomic states 
were initially populated: 



n st = |0 C | 2 < a + {t )a{t ) > +\<f) s \ 2 < b + (t )b(t ) > +2ite(## c < a(t )b(t ) >) (15) 

The functions <p c (t) and s (t) are c-f unctions that are independent solutions of the equa- 
tion (|I0"D. This equation with the initial conditions (O) appears also in the one mode case 
which was investigated in [JJ, (see appendix). Using the results of ]]J and substitute it into 
(|TT|) we find assimptotics for the self-consistent solution of our initial valued problem is as 
follows: 
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<f> s {t -> oo) = ^y / exp(2vrA) - i e ^-iS(NI)-<s(r)-i a r g r(iA) > 

c ( t _ oo) = e *A-KS(M)-iS(r) (16) 

where r = y/2j3t, S(t) = t 2 /A — Alnr and A = \ / y\ 2 /(2/3) Using this one can derive 

|0,(t-oc)| 2 = e 2 ^-l 
|0 c (t -> oo)| 2 = n sp = e 2?rA 

so if initially we have coherent atomic states \a > \(3 > where a\a >= a\a > and b\(3 >= 
P\P > then 



n st = \ a \ 2 e 2nX + \(3\ 2 {e 2nX - 1) + 2y/ e 2 ** - le* x \a\\/3\ x 
3n 

cos (_ + 2 S(\r \) + arg(r(*A)) + arg( 7 ) + arg(a) + arg(/3)) (18) 

so, depending on the initial arguments of a and (5 number of produced particles due to the 
stimulated transitions can be in the range n_ < n st < n + where 



n ± = {\a\e wX ±\(3\^/e 2 ^-l) 2 (19) 

Now, let us proceed to the problem of correlations between the two modes. We will 
investigate the case when initially there were no particles in atomic condensates. For this 
purpose we consider the quadrature phase amplitudes: 

X e (t) = (a{t) + b(t)y e + (a + (t) + b + (t))e^ e (20) 
If initially the atomic condensates are in vacuum states then: 

<X e \t)>=\Mt)e ie + ^s(t)e-' 9 \ 2 (21) 
and long time after curve crossing event we obtain: 

< X e \t -> OO) >= \y/e 2 * X - 1 + e ^- ! | L -*"g7-2.S(r)- l argr(iA)+2 l fl|2 ^ 

Choosing two orthogonal phase angles 9 + = \/2{^f + arg7 + 2S(r) + argT(iA)) and 
d_ = l/2(=z + arg 7 + 2S(t) + argT(iA)) we find: 



< X 2 ± >= \^/e 2 ^-l±e nX \ 2 (23) 



which is what one gets in the one mode case but with different interpretation. The 
result (|23|) means that dissociation forms atoms in a two modes squeezed state. On the 



applications of two mode squeezed states in atomic condensate see for example 0, |I5| , ||T6 

Finally we would like to discuss another application of our results. The same Hamiltonian 
(0) would appear in three components condensate (that differ for example by spin projection) 
if the following reaction during atomic collisions is possible: 



C + C — ► A + B (24) 

This process leads to a gipipa + b + term in the Hamiltonian which in the case of large number 
of C atoms can be substituted by g < ipip > a + b + . 

The time dependent chemical potential crossing can be achived by changing adiabatically 
the shape of the trap. This would produce the change of condensate concentration and hence 
of the chemichal potentials. 

If components A and B are distinct by some internal degree of freedom that is not 
conserved during the reaction then the Hamiltonian (|2|) should be generalized so that reaction 
go in all possible channels: 



H = Hi(t) = a + a + /i 2 (t)6 + 6 + 70a + ■y*a + b + + 7 a aa + 7*a + a + + ^bb + ^lb + b + (25) 

In this case even if we treat the field of initial condensate as a c-number, the resulting 
Landau-Zener problem in not analytically solvable, since all four operator equations become 
coupled. However, if additional symmetry exists, exact asymptotics can be found. For 
example, if the chemical potentials are equal due to some symmetry: 

ix x {t) = /x 2 (t) (26) 

then the equivalent Landau-Zener problem can be solved. For simplicity let us also assume 
that 7b = 7a = 7i- Then the operators evolution is govenered by the following equations: 

ia = ta + 76+ - 27^0+ 
ib = tb + 7*a + - 27*6+ 
this is accomplished by the hermitian conjugate equations. 
If we add the equations in (|27j) and denote c = a + b we get 
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(27) 



i'c = tc+ (7* - 2 7l *)c + (28) 

so that for c and c + we have a solvable two state Landau-Zener system. Accordingly, we 
have similar equation for operators d and d + , where d = a — b: 

id = td - (7* + 2jl)d + (29) 

At the end we can switch back to operators a and b. Such solution would be very unstable 
in respect of the terms that break the symmetry appear in the Hamiltonian given by ( |26| ) 
since there are states that have the same energy participating in the curve crossing. This 
instability was investigated in a set of articles fl7|l , [jTSfl , |L9 . 



In conclusion, we considered a time dependent problem of coherent dissociation of molec- 
ular condensate into distinct modes beyond the mean field approximation. We showed that it 
can be reduced to the case of two distinct atomic modes and molecular mode. We discussed 
in particular the case when the molecular condensate has zero momentum because in this 
case the interaction term with the two resulting modes is equal hence the coupling constant 
g between the molecular and the atomic field becomes 2g. This fact is very important since 
it leads to an enhancement in the number of the atomic condensate atoms that exist in the 
squeezed state. In the case of two modes condensate the corresponding multistate Landau- 
Zener problem can be solved exactly since the operators evolution equations decouple. The 
produced condensates are entangled in a two modes squeezed state. A more general case of 
the two mode condensate production during a curve crossing could also be exactly solvable 
if additional symmetry is imposed. However such system could be very unstable if terms 
that break the symmetry participate in the Hamiltonian of the system. 

APPENDIX 

The solutions of eq.(9) that satisfy the initial conditions (11) are given by 

W((j) 1} (j) 2 ) 

<Ps{t) = 2Z7 7777- (31) 

VK(0i,0 2 ) 

where W(<pi,<p2) = i{2(3 /^)exp(-Tx\/2) and 

8 



A = | T | 2 /(2/3) (32) 

is Landau-Zener parameter. The Functions <p\ and 02 are two standard solutions of parabolic 
cylinder equation with asymptotics at to — > — oo: 

Mto) ~ - exp(-^A + i- + iS^M)) (33) 
r 4 4 

2 (t o ) ~ -^exp (— jA - 4 " iS'(kol)), (34) 

and at £ — > +oo: 



Z I 7T 7T 

0i(*) ~ 1 — i \//3sinh (7rA) exp (—A — z iS{r) — iargTUX)) (35) 

|7| 4 2 

Mt) ~ -V^exp (^A - 4 - z5(|r |)) (36) 
7 ' 4 4 
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